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Abstract 

The lateral and longitudinal profiles of liadronic showers detected by a prototype 
of the ATLAS Iron-Scintillator Tile Hadron Calorimeter have been investigated. This 
calorimeter uses a unique longitudinal configuration of scintillator tiles. Using a fine- 
grained pion beam scan at 100 GeV, a detailed picture of transverse shower behavior 
is obtained. The underlying radial energy densities for four depth segments and for 
the entire calorimeter have been reconstructed. A three-dimensional hadronic shower 
parametrization has been developed. The results presented here are useful for under- 
standing the performance of iron-scintillator calorimeters, for developing fast simula- 
tions of hadronic showers, for many calorimetry problems requiring the integration of 
a shower energy deposition in a volume and for future calorimeter design. 

Keywords: Calorimetry; Computer data analysis. 
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1 Introduction 



We report on an experimental study of hadronic shower profiles detected by the prototype of 
the ATLAS Barrel Tile Hadron Calorimeter (Tile calorimeter) , . The innovative design 
of this calorimeter, using longitudinal segmentation of active and passive layers (see Fig. |I|), 
provides an interesting system for the measurement of hadronic shower profiles. Specifically, 
we have studied the transverse development of hadronic showers using 100 GeV pion beams 
and longitudinal development of hadronic showers using 20 - 300 GeV pion beams. 

Characteristics of shower development in hadron calorimeters have been published for 
some time. However, a complete quantitative description of transverse and longitudinal 
properties of hadronic showers does not exist [Q . The transverse profiles are usually expressed 
as a function of transverse coordinates, not the radius, and are integrated over the other 
coordinate @]. The three-dimensional parametrization of hadronic showers described here 
could be a useful starting point for fast simulations, which can be faster than full simulations 
at the microscopic level by several orders of magnitude ||^, 0, 0. 

The paper is organised as follows. In Section 2, the calorimeter and the test beam 
setup are briefiy described. In Section 3, the mathematical procedures for extracting the 
underling radial energy density of hadronic showers are developed. The obtained results on 
the transverse and longitudinal profiles, the radial energy densities and the radial contain- 
ment of hadronic shower are presented in Sections 4-7. Section 8 and 9 investigate the 
three-dimensional parametrization and electromagnetic fraction of hadronic shower. Finally 
Section 10 contains a summary and the conclusions. 

2 The Calorimeter 

The prototype Tile Calorimeter used for this study is composed of five modules stacked in 
the Y direction, as shown in Fig. 0. Each module spans 27r/64 in the azimuthal angle, 100 
cm in the Z direction, 180 cm in the X direction (about 9 interaction lengths. A/, or about 
80 effective radiation lengths, Xq), and has a front face of 100 x 20 cm^ The absorber 
structure of each module consists of 57 repeated "periods". Each period is 18 mm thick and 
consists of four layers. The first and third layers are formed by large trapezoidal steel plates 
(master plates) and span the full longitudinal dimension of the module. In the second and 
fourth layers, smaller trapezoidal steel plates (spacer plates) and scintillator tiles alternate 
along the X direction. These layers consist of 18 different trapezoids of steel or scintillator, 
each of 100 mm in depth. The master plates, spacer plates and scintillator tiles are 5 mm, 
4 mm and 3 mm thick, respectively. The iron to scintillator ratio is 4.67 : 1 by volume. The 
calorimeter thickness along the beam direction at the incidence angle (the angle between 
the incident particle direction and the normal to the calorimeter front face) of G = 10° 
corresponds to 1.49 m of iron equivalent [Q. 

Wavelength shifting fibers collect scintillation light from the tiles at both of their open 
(azimuthal) edges and bring it to photo-multipliers (PMTs) at the periphery of the calorime- 
ter. Each PMT views a specific group of tiles through the corresponding bundle of fibers. 
The modules are divided into five segments along Z. They are also longitudinally segmented 
(along X) into four depth segments. The readout cells have a lateral dimensions of 200 
mm along Z, and longitudinal dimensions of 300, 400, 500, 600 mm for depth segments 1 
- 4, corresponding to 1.5, 2, 2.5 and 3 A/ at 9 = 0° respectively. Along Y, the cell sizes 
vary between about 200 and 370 mm depending on the X (depth) coordinate. We recorded 
energies for 100 different cells for each event P). 

The calorimeter was placed on a scanning table that allowed movement in any direction. 
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Upstream of the calorimeter, a trigger counter telescope (SI - S3) was installed, defining a 
beam spot approximately 20 mm in diameter. Two delay- line wire chambers (BCl - BC2), 
each with [Z, Y) readout, allowed the impact point of beam particles on the calorimeter face 



to be reconstructed to better than ±1 mm |10]. "Muon walls" were placed behind (800x800 



mm^) and on the positive Z side (400x1150 mm^) of the calorimeter modules to measure 
longitudinal and lateral hadronic shower leakage ||11||. 



We used the TILEMON program |0 to convert the raw calorimeter data into PAW 
Ntuples |l^ containing calibrated cell energies and other information used in this study. 

The data used for the study of lateral profiles were collected in 1995 during a special 
Z-scan run at the CERN SPS test beam. The calorimeter was exposed to 100 GeV negative 
pions at a 10° angle with varying impact points in the Z-range from —360 to +200 mm. A 
total of > 300,000 events have been analysed; for the lateral profile study only events without 
lateral leakage were used. The uniformity of the calorimeter's response for this Z-scan is 



estimated to be 1% 114 



The data used for the study of longitudinal profiles were obtained using 20 - 300 GeV 
negative pions at a 20° angle and were also taken in 1995 during the same test beam run. 



3 Extracting the Underlying Radial Energy Density 

In this investigation we use a coordinate system based on the incident particle direction. 
The impact point of the incident particle at the calorimeter front face defines the origin of 
coordinate system. The incident particle direction forms the x axis, while the y axis is in the 
same direction as Y defined in Section 2. The normal to the xy surface defines the z axis. 

We measure the energy deposition in each calorimeter cell for every event. In the ijk-cel\ 
of the calorimeter with the volume Vijk and cell center coordinates {xc,yc,Zc), the energy 
deposition Eij^ is 

^ ...... ,1) 



V,, 



ijk 



where f(x,y,z) is the three-dimensional hadronic shower energy density function. Due to 
the azimuthal symmetry of shower profiles, the density f{x,y,z) is only a function of the 
radius r = \^y^~+^ from the shower axis and the longitudinal coordinate x. Then 

Eijk{xc,yc,Zc) = J ^(x,r) rdrdcpdx, (2) 



ijk 



where is the azimuthal angle and "^{x, r) has the form of a joint probability density function 
(p.d.f.) [|T^. The joint p.d.f. can be further decomposed as a product of the marginal p.d.f., 
dE{x)/dx, and the conditional p.d.f., <l>(x, r), 

\ij{x,r) = '^^j^^ ■ ^{x,r). (3) 

(JjJb 

The longitudinal density dE/dx is defined as 

dE{x) 



dx 



j j f{x, y, z) dydz. (4) 



Finally, the radial density function $(r) for a given depth segment is 

$(r)=Eo<l'(a;,r), (5) 
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where Eq is the total shower energy deposition into the depth segment for fixed x. 

There are several methods for extracting the radial density $(r) from the measured 
distributions of energy depositions. One method is to unfold $(r) using expression (|^). This 
method was used in the analysis of the data from the lead-scintillating fiber calorimeter 
ll6| . Several analytic forms of $(r) were tried, but the simplest that describes the energy 
deposition in cells was a combination of an exponential and a Gaussian: 

^(r) = — e '^1 + — e , (6) 

where bi and /ij are the free parameters. 

Another method for extracting the radial density is to use the marginal density function 

OO X2 

f\z) = j j fix,y,z) dxdy. (7) 

— OO Xl 

which is related to the radial density $(r) 



OO 



^ 2 / ^^t^L. (8) 



This method was used for extracting the electron shower transverse profile from the 
GAMS-2000 electromagnetic calorimeter data |18|. The above integral equation (§) can be 
reduced to an Abelian equation by replacing variables In the following solution to 
equation was obtained 

^/ \ _ I d J f{z) dz^ 

^''^ ^ ~n d^ J ^ ^ 

9 

For our study, we used the sum of three exponential functions to parameterize f{z) as 

= ^ E «^ e"^' (10) 
i=i 

where z is the transverse coordinate, Eq, a^, Aj are free parameters, B = J2l=i cuK, J2l=i = 
1 and f{z)dz = Eq. In this case, the radial density function, obtained by integration 
and differentiation of equation (^, is 

where Kq is the modified Bessel function. This function goes to oo as r — > and goes to 
zero as r ^ OO. 

We define a column of five cells in a depth segment as a tower. Using the parametrization 
shown in equation (0), we can show that the energy deposition in a tower [^, E{z) = 
Iz-h/2 fi^)(^^^ can be written as 



E{z) = ^0 - f E a,A, cosh(g) e for \z\ < ^, (12) 

i=l ^ I 

3 1^1 

E{z) = f E «^Ai sinh(2^) e"^, for \z\ > -, (13) 

i=l 



2' 
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where h is the size of the front face of the tower along the z axis. Note that as /i — 0, we 
get E{z)/h f{z). As /i ^ oo, we find that E{0) Eq. 

The full width at half maximum (FWHM) of an energy deposition profile for small values 
of {FWHM ~ h)/h can be approximated by 

FWHM = h+ 2E(V2)-E(0) i^^j 
f E Oi sinh(5|-) e 

1=1 ' 

We will show below that this approximation agrees well with our data. 
A cumulative function may be derived from the density function as 

z 

F{z) = J f{z) dz. (15) 

— oo 

For our parametrization in equation (p^, the cumulative function becomes 

F{z) = E a,Xi et, for z < 0, (16) 

i=l 

F{z)= ^o-§ E aiXie'^, farz>0, (17) 

i = l 

where z is the position of the edge of a tower along the z axis. Note that the cumulative 
function does not depend on the cell size h. We can construct the cumulative function and 
deconvolute the density f(z) from it for any size calorimeter cell. Note also that cumulative 
function is well behaved at the key points: F(— oo) = 0, F{0) = Eq/2, and F(oo) = Eq. 
The radial containment of a shower as a function of r is given by 

/(r) = J J $(r) rdrd(j) = - ^ a, i^i (^^ 

' 

where Ki is the modified Bessel function. As r —>■ oo, the function rKi{r) tends to zero and 
we get /(oo) = £"0, as expected. 

We use two methods to extract the radial density function $(r). One method is to unfold 
$(r) from (^. Another method is to use the expression (|^) after we have obtained the 
marginal density function f{z). There are three ways to extract f{z): by fitting the energy 
deposition E{z) [^], by fitting the cumulative function F{z) and by directly extracting f{z) 
by numerical differentiation of the cumulative function. The effectiveness of these various 
methods depend on the scope and quality of the experimental data. 



4 Transverse Behaviour of Hadronic Showers 

Figure ^ shows the energy depositions in towers for depth segments 1 - 4 as a function of the 
z coordinate of the center of the tower. Figure ^ shows the same for the entire calorimeter 
(the sum of the histograms presented in Fig. 

These Figures are profile histograms and give the energy deposited in any tower for 
all analysed events in bins of the z coordinate. Here the coordinate system is linked to the 
incident particle direction where ^ = is the coordinate of the beam impact points at the 
calorimeter front face. Figure |^ schematically shows a top view of the experimental setup 
and indicates the z coordinate of each tower. The zi, Z2, z^ and Z4 are the distances between 
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the centre of towers (for the four depth segments) and the direction of beam particle. The 
values of the 2;-coordinate of the tower centers are negative to the left of the beam, positive 
to the right of the beam and range from —750 mm to +600 mm. To avoid edge effects, we 
present tower energy depositions in the range from —650 mm to +500 mm. Note that the 
total tower height (about 1.0 m at the front face, and about 1.8 m at the back) is sufficient 
for shower measurements without significant leakage in the vertical direction. 

As mentioned earlier, events with significant lateral leakage (identified by a clear minimum- 
ionising signal in the lateral muon wall) were discarded. The resulting left-right asymmetries 
in the distributions of Figures ^ and |^ are very small. 

As will be shown later (Section 6), the 99% containment radius is less than 500 mm. 

The fine-grained 2;-scan provided many different beam impact locations within the ca- 
lorimeter. Due to this, we obtained a detailed picture of the transverse shower behaviour 
in the calorimeter. The tower energy depositions shown in Figures ^ and ^ span a range 
of about three orders of magnitude. The plateau for \z\ < 100 mm {h/2) and the fall-off 
at large \z\ are apparent. Similar behaviour of the transverse profiles was observed in other 



calorimeters as well [161, 21 



We used the distributions in Figs. |^ and ^ to extract the underlying marginal densities 
function for four depth segments of the calorimeter and for the entire calorimeter. The 
solid curves in these figures are the results of the fit with equations (12) and (^). The fits 
typically differ from the experimental distribution by less than 5%. 

In comparison with |^ and |2^ , where the transverse profiles exist only for distances less 
than 250 mm, our more extended profiles (up to 650 mm) require that the third exponential 
term be introduced. The parameters Oj and Aj, obtained by fitting, are listed in Table 0. The 
values of the parameter i?o, the average energy shower deposition in a given depth segment, 
are listed in Table ^. 

We have compared our values of Ai and A2 with the ones from the conventional iron- 



scintillator calorimeter described in |2^. At 100 GeV, our results for the entire calorimeter 



are 23 + 1 mm and 58 + 4 mm for Ai and A2 respectively. They agree well with the ones 
from [Q, which are 18 + 3 mm and 57 + 4 mm. 

We determined the FWHM of energy deposition profiles (Figs. § and 0) using formula 
(|T^). The characteristic FWHM are found to be approximately equal to transverse tower 
size. The relative difference of FWHM from transverse tower size, {FWHM — h)/h, amount 
to 2% for depth segment 1, depth segment 2 and for the entire calorimeter, 7% for depth 
segment 3 and 15% for depth segment 4. 

Figure |^ shows the calculated marginal density function f{z) and the energy deposition 
function, E{z)/h, at various transverse sizes of tower h = 50, 200, 300 and 800 mm using 
the obtained parameters for the entire calorimeter. As a result of the volume integration, the 
sharp f{z) is transformed to the wide function E[z)/h, which clearly shows its relationship 
to the transverse width of a tower. The values of FWHM are 40 mm for f{z) and 204 mm 
for E{z)/h a.t h = 200 mm. Note that the transverse dimensions of a tower vary from 300 
mm to 800 mm for the different depth segments in final ATLAS Tile calorimeter design. The 
difference between /(O) and E{0)/h becomes less then 5% only at h less then 6 mm. 

The parameters and Ai as a function of x (in units of Af^) are displayed in Figs. ^ 
and ^ Here A^*^ = 207 mm is the nuclear interaction length for pions in iron In these 
calculations, the effect of the 10° incidence beam angle has been corrected. As can be seen 
from Figs. |^ and |^, the value of ai decreases and the values of the remaining parameters, 
02,03 and Aj, increase as the shower develops. This is a reflection of the fact that as the 
hadronic shower propagates into the calorimeter it becomes broader. Note also that the a, 
and Aj parameters demonstrate linear behaviour as a function of x. The lines shown are flts 
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to the linear equations 

ai{x) = ctj + PiX (19) 

and 

Ai(x) = 7i + 5iX. (20) 

The values of the parameters f3i, 7^ and 6i are presented in Table |^. It is interesting to 
note that the linear behaviour of the slope exponential was also observed for the low-density 
fine-grained flash chamber calorimeter However, some non-linear behaviour of the slope 
of the halo component was demonstrated for the uranium-scintillator ZEUS calorimeter at 
interaction lengths more than 5 A/ P3 . 



Similar results were obtained for the cumulative function distributions. The cumulative 
function F{z) is given by 

F{z) = j:F\z), (21) 

k=l 

where F^{z) is the cumulative function for depth segment k. For each event, F^{z) is 



''•max 
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F\z)=Y.T.E^,k, (22) 
i=i j=i 

where imax = 1, . • • , 5 is the last tower number in the sum. 

Figures and [TO present the cumulative functions for four depth segments and for the 



entire calorimeter. The curves are fits of equations (16) and (|T^) to the data. Systematic 
and statistical errors are again added in quadrature. The results of the cumulative function 
fits are less reliable and in what follows we use the results from energy depositions in a tower. 

However, the marginal density functions determined by the two methods (by using the 
energy deposition spectrum and the cumulative function) are in reasonable agreement. 



5 Radial Hadronic Shower Energy Density 

Using formula (|ll]) and the values of the parameters Oj, A^, given in Table |l|, we have 
determined the underlying radial hadronic shower energy density functions, $(r). The results 
are shown in Figure for depth segments 1-4 and in Figure ^ for the entire calorimeter. 
The contributions of the three terms of $(r) are also shown. 

The functions $(r) for separate depth segments of a calorimeter are given in this paper for 
the first time. The function $(r) for the entire calorimeter was previously given for a Lead- 
scintillating fiber calorimeter it is given here for the first time for an Iron-scintillator 
calorimeter. The analytical functions giving the radial energy density for different depth 
segments allow to easily obtain the shower energy deposition in any calorimeter cell, shower 
containment fractions and the cylinder radii for any given shower containment fraction. 

The function $(r) for the entire calorimeter has been compared with the one for the lead- 



scintillating fiber calorimeter of ref. |16], that has about the same effective nuclear interaction 



length for pions (namely 251 mm for the tile and 244 mm for the fiber calorimeter [|20 



The two radial density functions are rather similar as seen in Fig. [T^. The lead-scintillating 
fiber calorimeter density function $(r), which was obtained from a 80 GeV n~ grid scan at 
an angle of 2° with respect to the fiber direction, was parametrized using formula (|^) with 
bi = 0.169 pC/mm, 62 = 0.677 pC/mm, fii = 140 mm and fi2 = 42.4 mm. For the sake of 
comparing the radial density functions of the two calorimeters, the distribution from JIB 



was 



normalised to the $(r) of the Tile calorimeter. Precise agreement between these functions 
should not be expected because of the effect of the different absorber materials used in the 
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two detectors (e. g. the radiation/interaction length ratio for the Tile calorimeter is three 
times larger than for lead-scintillating fiber calorimeter [^), the values of e/h are different, 



as is hadronic activity of showers because fewer neutrons are produced in iron than in lead 



6 Radial Containment 

An other issue on which new results are presented here is the longitudinal development of 
shower transverse dimensions. The parametrization of the radial density function, $(r), 
was integrated to yield the shower containment as a function of the radius, I{r). Figure p!4| 
shows the transverse containment of the pion shower, /( , clS db function of r for four depth 
segments and for the entire calorimeter. 

In Table |] and Fig. 15 the radii of cylinders for the given shower containment (90%, 
95% and 99%) extracted from Fig. [1^ as a function of depth are shown. The centers of 
depth segments, x, are given in units of X^^. Solid lines are the linear fits to the data: 
r(90%) = (85±6) + (37±3)a;, r(95%) = (134±9) + (45±3)x, r(99%) = (349±7) + (22±2)a; 
(mm). As can be seen, these containment radii increase linearly with depth. Such a linear 
increase of 95% lateral shower containment with depth is also observed in an other iron- 



scintillator calorimeter at 50 and 140 GeV |26|. It is interesting to note that the shower 



radius for 95% radial containment for the entire calorimeter is equal to X"^^^ = 251 mm 
which justifies the frequently encountered statement that r(95%) ~ Xj [Q, where A/ is X^^^^ 
in our case. For the entire Tile calorimeter the 99% containment radius is equal to 1.7 ± 0.1 

Based on our study, we believe that it is a poor approximation to regard the values ob- 
tained from the marginal density function or the energy depositions in strips as the measure 
of the transverse shower containment, as was done in In that paper the value of 1.1 A^-^-^ 
was obtained for 99% containment at 100 GeV, and the conclusion was drawn that their 
"result is consistent with the rule of thumb that a shower is contained within a cylinder of 
radius equal to the interaction length of a calorimeter material" . However Tile calorimeter 
measurements show that the cylinder radius for 99% shower containment is about two in- 
teraction lengths. If we extract the lateral shower containment dimension using instead the 
integrated function F{z), given in Fig. |10|, we obtain the value of 300 mm or 1.2 X"^^-^ , which 
agrees with [Q. 

7 Longitudinal Profile 

We have examined the differential deposition of energy AE/Ax as a function of x, the 
distance along the shower axis. Table ^ lists the centers in x of the depth segments, x, 
and the lengths along x of the depth segments. Ax, in units of A^^, the average shower 
energy depositions in various depth segments, Eq, and the energy depositions per interaction 
length A^^, AE/Ax. Note that the values of Eq have been obtained taking into account the 



longitudinal energy leakage which amounts to 1.8 GeV for 100 GeV |n 



Our values of AE/Ax together with the data of |2^ and Monte Carlo predictions 



(GEANT-FLUKA + MICAP) || are shown in Fig. [T|. The longitudinal energy depo- 
sition for our calorimeter using longitudinal orientation of the scintillating tiles is in good 
agreement with that of a conventional iron-scintillator calorimeter. 

The longitudinal profile, AE/Ax, may be approximated using two parametrizations. The 
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first form is 



dE{x) _ Ef f3 



a+l 



dx 



r(a + 11 



a;" e 



-I3x 



(23) 



where Ef = Ebeam, and a and are free parameters. Our data at 100 GeV and those of Ref. 



27| at 100 GeV were jointly fit to this expression; the fit is shown in Fig. [16 



The second form is the analytical representation of the longitudinal shower profile from 
the front of the calorimeter 

dE{x) 



+ 



(l-w) \j 



X 



X 



iFi(l, a + l, (^^-1) ^ 



(24) 



where iFi is the confluent hypergeometric function |g9|. Here the depth variable, x, is the 
depth in equivalent Fe, Xq is the radiation length in Fe and in this case A/ is A^*^. The 
normalisation factor, N, is given by 



N = E, 



beam 



dE{x] 
dx 



dx 



E, 



beam 



Xi r(a) {w Xo b-^ + {l-w) A/ rf-«) ' 



(25) 



This form was suggested in and derived by integration over the shower vertex positions 
of the longitudinal shower development from the shower origin 



dE{x) 
dx 



dEJx 



dx 



(26) 



where x^ is the coordinate of the shower vertex. (This is necessary because with the Tile 
calorimeter longitudinal segmentation the shower vertex is not measured). For the parame- 
trization of longitudinal shower development, the well known parametrization suggested by 
Bock et al. |^ has been used 



dEsjx) 
dx 



N Iw — 



X 



a-l 



-fe— 

^0 



+ (1 



W 



X 

a7 



a-l 



(27) 



where a, b, d, w are parameters. 

We compare the form ( ^If ) to the experimental points at 100 GeV using the parameters 
calculated in Refs. ||^ and [^. Note that now we are not performing a fit but checking 



how well the general form (^4]) together with two sets of parameters for iron-scintillator 
calorimeters describe our data. As shown in Fig. [TO, both sets of parameters work rather 
well in describing the 100 GeV data. 

Turning next to the longitudinal shower development at different energies, in Fig. ^ our 
values of AE/Ax for 20 - 300 GeV together with the data from [2? are shown. The solid 
and dashed lines are calculations with function ( p4|) using parameters from and re- 



spectively. Again, we observe reasonable agreement between our data and the corresponding 
data for conventional iron-scintillator calorimeter on one hand, and between data and the 
parametrizations described above. Note that the fit in has been performed in the energy 
range from 10 to 140 GeV; hence the curves for 200 and 300 GeV should be considered as 
extrapolations. It is not too surprising that at these energies the agreement is significantly 
worse, particularly at 300 GeV. In contrast, the parameters of |^ were derived from data 
spanning the range 15 - 400 GeV, and are in much closer agreement with our data. 
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8 The parametrization of Hadronic Showers 



The three-dimensional parametrization for spatial hadronic shower development is 

3 I \ 

V Q»W TV- / r \ 



2tt J2 ai{x)Xi{x) 



i=l 



where dE{x)/dx, defined by equation (0), is the longitudinal energy deposition, the func- 
tions ai{x) and Aj(x) are given by equations (pJ]) and (pOD, and Kq is the modified Bessel 
function. 

This explicit three dimensional parametrization can be used as a convenient tool for many 
calorimetry problems requiring the integration of a shower energy deposition in a volume 
and the reconstruction of the shower coordinates. 



9 Electromagnetic Fraction of Hadronic Showers 

One of the important issues in the understanding of hadronic showers is the electromagnetic 
component of the shower, i. e. the fraction of energy going into 7r° production and its depen- 
dence on radial and longitudinal coordinates, /^o(r, x). Following [|1^], we assume that the 
electromagnetic part of a hadronic shower is the prominent central core, which in our case is 
the first term in the expression ([TI| ) for the radial energy density function, $(r). Integrating 
/ttO over r we get 

/.o = (29) 

i=l 

For the entire Tile calorimeter this value is (53 ±3)% at 100 GeV. 

The observed tt" fraction, /^o, is related to the intrinsic actual fraction, f^o, by the 
equation 

^^^^ ~eEi^ + hE',- {e/h - 1) • /:o(E) + 1' ^'"^ 

where E'^^ and i?^ are the intrinsic electromagnetic and hadronic parts of shower energy, 
e and h are the coefficients of conversion of intrinsic electromagnetic and hadronic energies 
into observable signals, /^o = £'em/(-^em + ^h)- 

There are two analytic forms for the intrinsic 7r° fraction suggested by Groom 

flo{E) = l-[^^j (31) 



and Wigmans 

where E'q = 1 GeV, m = 0.85 and k = 0.11. We calculated /^o using the value e/h 



fl,[E) = k.ln[^Y (32) 



1.34 it 0.03 for our calorimeter p[, and obtained the curves shown in Fig. 

Our result at 100 GeV is compared in Fig. |1^ to the modified Groom and Wigmans 
parametrizations and to results from the Monte Carlo codes CALOR [^, GEANT-GEISHA 
^ and GEANT-CALOR (the latter code is an implementation of CALOR89 differing 
from GEANT-FLUKA only for hadronic interactions below 10 GeV). Note that the Monte 
Carlo calculations were performed for the intrinsic 7r° fraction, f^o{E), and therefore the 
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results were modified by us according to (pOl). As can be seen from Fig. |18|, our calculated 
value of /^o is about one standard deviation lower than two of the Monte Carlo results and 
the Groom and Wigmans parametrizations. 

Figure 19 shows the fractions fnoir) as a function of r. As can be seen, the fractions 
f-Ko{r) for the entire calorimeter and for depth segments 1-3 amount to about 90% as r — >• 
and decrease to about 1% as r ^ Kl^ ■ However for depth segment 4 the value of f-Ko{r) 
amounts to only 50% as r — > and decreases slowly to about 10% as r — X^J^ ■ 

Figure ^ shows the values of f-,^o{x) as a function of a;, as well as the linear fit which 
gives Uo{x) = (75 ± 2) - (8.4 ± 0.4)a; (%). 

Using the values of f-^o (x) and energy depositions for various depth segments, we obtained 
the contributions from the electromagnetic and hadronic parts of hadronic showers in Fig. 



16] . The curves represent a fit to the electromagnetic and hadronic components of the shower 



using equation (pS]). Ej is set equal to f-^oEheam for the electromagnetic fraction and (1 — 
/^o)£^feeam for the hadronic fraction. The electromagnetic component of a hadronic shower 
rise and decrease more rapidly than the hadronic one (ogm = 1.4 ±0.1, ah = 1.1 ±0.1, /3em = 
1.12 ±0.04, Ph = 0.65 ±0.05). The shower maximum position {xmax = (o^/P) Kl^) occurs at 
a shorter distance from the calorimeter front face {x^rnax = 1-23 Kl^ ^ ^max = 1-85 A^-^-^). At 
depth segments greater than 4 X'^J^ ., the hadronic fraction of the shower begins to dominate. 
This is natural since the energy of the secondary hadrons is too low to permit significant 
pion production. 



10 Summary and Conclusions 

We have investigated the lateral development of hadronic showers using 100 GeV pion beam 
data at an incidence angle of O = 10° for impact points z in the range from —360 to 200 
mm and the longitudinal development of hadronic showers using 20 - 300 GeV pion beams 
at an incidence angle of G = 20°. 

Some useful formulae for the investigation of lateral profiles have been derived using a 
three-exponential form of the marginal density function f{z). 

We have obtained for four depth segments and for the entire calorimeter: energy deposi- 
tions in towers, E{z)\ cumulative functions, F{z)\ underlying radial energy densities, $(r); 
the contained fraction of a shower as a function of radius, /(r); the radii of cylinders for a 
given shower containment fraction; the fractions of the electromagnetic and hadronic parts 
of a shower; differential longitudinal energy deposition AE / ^x] and a three-dimensional 
hadronic shower parametrization. 

We have compared our data with those from a conventional iron-scintillator calorimeter, 
those from a lead-scintillator fiber calorimeter, and with Monte Carlo calculations. We have 
found that there is general reasonable agreement in the behaviour of the Tile calorimeter 
radial density functions and those of the lead-scintillating fiber calorimeter; that the longi- 
tudinal profile agrees with that of a conventional iron-scintillator calorimeter and the Monte 
Carlo predictions; that the value at 100 GeV of the calculated fraction of energy going into 7r° 
production in a hadronic shower, 7,^0, agrees with the Groom and Wigmans parametrizations 
and with some of the Monte Carlo predictions. 

The three-dimensional parametrization of hadronic showers that we obtained allows direct 
use in any application that requires volume integration of shower energy depositions and 
position reconstruction. The experimental data on the transverse and longitudinal profiles, 
the radial energy densities and the three-dimensional hadronic shower parametrization are 
useful for understanding the performance of the Tile calorimeter, but might find broader 
application in Monte Carlo modeling of hadronic showers, in particular in fast simulations. 
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and for future calorimeter design. 
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Table 1: The parameters and Aj obtained by fitting the transverse shower profiles for 
four depth segments and the entire calorimeter at 100 GeV. 



Depth 






Ai (mm) 


02 


A2 (mm) 


03 


A3 (mm) 


1 


0.6 


0.88 ±0.07 


17± 2 


0.12 ±0.07 


48 ± 14 


0.004 ±0.002 


430 ± 240 


2 


2.0 


0.79 ±0.06 


25 ±2 


0.20 ±0.06 


52 ±6 


0.014 ±0.006 


220 ± 40 


3 


3.8 


0.69 ±0.03 


32 ±8 


0.28 ±0.03 


71 ± 13 


0.029 ±0.005 


280 ± 30 


4 


6.0 


0.41 ±0.05 


51 ±10 


0.52 ±0.06 


73± 18 


0.07 ±0.03 


380 ± 140 


all four 




0.78 ± 0.08 


23 ± 1 


0.20 ±0.08 


58 ±4 


0.015 ± 0.004 


290 ± 40 



Table 2: The values of the parameters ai, f3i, 7^ and 5i 





ai 


A (i/A.) 




7i (mm) 


5i (mm/A^) 


ai 


0.99 ±0.06 


-0.088 ±0.015 


Ai 


13 ±2 


6± 1 


a2 


0.04 ± 0.06 


0.071 ± 0.015 


A2 


42± 10 


6±4 


as 


-0.001 ± 0.002 


0.008 ± 0.002 


A3 


170 ± 80 


29 ±23 
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Table 3: The radii of cylinders for the given shower containment at 100 GeV. 







90% 


95% 


99% 


0.6 


0.44 


0.64 


1.43 


2.0 


0.60 


0.88 


1.55 


3.8 


0.92 


1.27 


1.75 


6.0 


1.24 


1.55 


1.87 


all four 


0.72 


1.04 


1.67 



Table 4: Average energy shower depositions at various depth segments at 100 GeV. 







Eo (GeV) 


AE/Ax (GeV/A^'^) 


0.6 


1.2 


25.0±0.3 


20.8±0.3 


2.0 


1.6 


42.8±0.2 


26.8±0.1 


3.8 


2.0 


22.0±0.1 


ll.OiO.l 


6.0 


2.4 


8.4±0.5 


3.4±0.2 
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Muonwall 5 Module Detector 

Not to scale 



Figure 2: Schematic layout of the experimental setup. SI - S3 are beam trigger scintillators, 
and BCl - BC2 are (Z,Y) proportional chambers. 
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Figure 3: Energy depositions of 100 GeV pions in towers of depth segments 1 - 4 as a 
function of the z coordinate: top left is for depth segment 1, top right is for depth segment 
2, bottom left is for depth segment 3, bottom right is for depth segment 4. Only statistical 
errors are shown. Curves are fits of equations (12) and (0) to the data. 
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Figure 4: Energy depositions in towers, summed over all calorimeter depth segments, as a 
function of the z coordinate. Only statistical errors are shown. The curve is the result of 



the fit by formulas (12) and (13 
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Figure 5: Schematic layout (top view) of Tile calorimeter experimental setup, zi - zi^ are 
the distances between the centre of towers (for the four depth segments) and the direction 
of the beam particle. 
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Figure 6: The calculated marginal density function f{z) (the solid line) and the energy depo- 
sition function, E{z)/h, for various transverse sizes of a tower (h) : 50 mm (the dash-dotted 
line), 200 mm (the dashed line), 300 mm (the thick dotted line), 800 mm (the thin dotted 
line). The parameters for the entire calorimeter (see Table |1|) are used in the calculations. 
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Figure 7: X dependences of the parameters Oj: the triangles are the ai parameter, the 
diamonds are the 02, the squares are the 03. 
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Figure 8: X dependences of the parameters Aj: the triangles are the Ai parameter, the 
diamonds are the A2, the squares are the A3. 
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Figure 9: Cumulative functions F{z) for depth segments 1 - 4 as a function of the z 
coordinate: top left is for depth segment 1, top right is for depth segment 2, bottom left is 
for depth segment 3, bottom right is for depth segment 4. Statistical and systematic errors, 
summed in quadrature, are shown. Curves are fits of equations (16) and (|1^) to the data. 
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Figure 10: The cumulative function F{z) for the entire calorimeter as a function of the z 
coordinate. Only statistical errors are shown. Curves are fits of equations (16) and ([TTp to 
the data. 
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Figure 11: Radial energy density, $(r), as a function of r for Tile calorimeter for depth 
segments 1-4: top left is for depth segment 1, top right is for depth segment 2, bottom left 
is for depth segment 3, bottom right is for depth segment 4. The solid lines are the energy 
densities $(r), the dashed lines are the contribution from the first term from equation ([TTD , 
the dash-dotted lines are the contribution from the second term, the dotted lines are the 
contribution from the third term. 



27 




Figure 12: The radial energy density as a function of r (in units of A^'^'^) for Tile calorimeter 
(the solid line), the contribution to $(r) from the first term in equation ([TI|) (the dashed 
line), the contribution to $(r) from the second term (the dash-dotted line), the contribution 
to $(r) from the third term (the dotted line). 
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Figure 13: Comparison of the radial energy densities as a function of r (in units of \%^^) 
for Tile calorimeter (the solid hne) and lead-scintillating fiber calorimeter (the dash-dotted 
line) . 
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Figure 14: Containment of shower /(r) (the sohd hne) as a function of radius for the entire 
Tile calorimeter. The dash-dotted hne is the contribution from the first depth segment, the 
dashed line is the contribution from the second depth segment, the thin dotted line is the 
contribution from the third depth segment, the thick dotted line is the contribution from 
the fourth depth segment. 
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Figure 16: The longitudinal profile (circles) of the hadronic shower at 100 GeV as a function 
of the longitudinal coordinate x in units of A^^. Open triangles are data from the calorimeter 
of Ref. diamonds are the Monte Carlo (GEANT-FLUKA) predictions. The dash-dotted 
line is the fit by function (|23|) , the solid line is calculated with function (|2^) with parameters 
from Ref. [^|, the dashed line is calculated with function (24) with parameters from Ref. f^. 
The electromagnetic and hadronic components of the shower (crosses and squares), together 
with their fits using (|23|) , are discussed in Section 9. 
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Figure 17: Longitudinal profiles of the hadronic showers from 20 GeV (open stars), 50 GeV 
(open squares) and 100 GeV (open circles) pions as a function of the longitudinal coordinate 
X in units of A/ for a conventional iron-scintillator calorimeter |^ and of 20 GeV (black 
stars), 50 GeV (black squares), 100 GeV (black circles), 150 GeV (black up triangles), 200 
GeV (asterisks), 300 GeV (black down triangles) for pions at 20° and of 100 GeV (black 
circles) for pions at 10° for Tile iron-scintillator calorimeter. The solid lines are calculated 



with function (P4 ) with parameters from ||2^. The dashed lines are calculated with function 
with parameters from M. 
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Figure 18: The fraction /^o in hadronic showers versus the beam energy. The star is our data, 
the sohd curve is the Groom parametrization, the dashed curve is the Wigmans parametriza- 
tion, squares are the GEANT-CALOR predictions, circles are the GEANT-GHEISHA pre- 
dictions and crosses are the CALOR predictions. 
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Figure 19: The f^o fractions of hadronic showers as a function of radius. The sohd hue 
is the fT^o{r) for the entire Tile calorimeter. The dash-dotted hne is the f.,^o(r) for the first 
depth segment, the dashed line is the fno{r) for the second depth segment, the thin dotted 
line is the /,r"(^") for the third depth segment, the thick dotted line is the /7ro(r) for the fourth 
depth segment. 
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Figure 20: The /7ro(x) fractions of hadronic showers as a function of x. 
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